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THE  IMPLEMENTATION  OF  THE  PLOTTING  PROGRAM  PLOTEZ 


1 .  INTRODUCTION 


PLOTEZ  is  a  plotting  program  which  produces  graphical  output  of  data 
from  a  VAX  computer.  It  utilizes  the  Tektronix  PLOT1 0  IGL  library  of 
graphics  subroutines  [1]  to  create  two-dimensional  plots  of  data  on  Tektronix- 
type  graphics  terminals.  The  horizontal  (X)  axis  usually  represents  the 
independent  data  (X-values)  while  the  vertical  (Y)  axis  represents  the 
dependent  (Y-values)  data  set.  It  should  be  emphasized  that  PLOTEZ  is  a 
plotting  program  and  not  a  graphics  package;  that  is,  no  manipulation  of  the 
data  plotted  on  the  screen  is  possible. 


PLOTEZ  has  had  a  rich  history.  it  had  its  origins  as  the  simple 
plotting  program  XYPLOT  used  for  plotting  PHA  (Pulse  Height  Analysis)  gamma- 
ray  spectra  on  a  PDP-15  minicomputer  located  in  the  Physics  School  at 
Melbourne  University  in  the  early  1970's.  The  output  device  used  was  a 
Tektronix  611  storage  display  scope  and  graphical  output  achieved  using  some 
primitive  pen-up,  pen-down  commands  which  came  supplied  with  the  PDP-15 
system.  Hard  copy  was  achieved  by  doing  a  screen-dump  using  a  photocopy- type 
device. 


With  the  introduction  of  a  VAX-11/780  computer  in  1979  at  Melbourne 
University  this  program  was  ported  to  the  VAX,  considerably  enhanced  (and 
renamed  PLOT).  Graphical  output  to  TEK-4010  terminals  was  achieved  by 
sending  special  character  strings  to  the  terminal  to  put  it  into  graphics 
mode,  whereas  text  output  (for  numbers  on  the  axes  and  titles)  was  achieved  by 
simply  exiting  from  graphics  mode  and  using  formatted  FORTRAN  WRITE 
instructions.  This  was  the  only  way  to  proceed  at  the  time  since  the  VAX  had 
no  plotting  package  installed.  PLOT  ran  in  this  fashion  for  several  years  on 
the  VAX-VMS  system  very  successfully  and  has  now  been  adapted  to  run  on  the 
Data  General  MV-8000  in  the  Physics  School  at  Melbourne  University. 


In  1981 ,  the  author  accepted  a  position  at  the  University  of 
Massachussetts  at  Amherst  (USA)  and  found  that  there  was  no  suitable  plotting 
program  for  the  RSX-1 1  system  running  on  the  PDP-11/35  which  the  Nuclear 


Physics  group  used  for  data  analysis.  The  PLOT  program  was  then  converted  to 
run  on  this  machine  using  the  early  PL0T1 0  Terminal  Control  System  (TCS) 
graphics  library  [2].  Later,  in  1983,  when  the  group  upgraded  their  computer 
facility  to  a  VAX-1 1/750  the  code  was  adapted  to  run  on  this  machine.  This 
program  was  very  successful  and  used  widely.  It  is  still  running  at  the 
University  of  Massachussetts  and  elsewhere.  Among  the  modifications  applied 
at  UMass  was  the  ability  to  produce  semi-log  plots  which  were  needed  for 
plotting  electron  scattering  cross  sections.  These  typically  varied  by  orders 
of  magnitude  for  small  changes  in  momentum  transferred  to  the  nucleus. 


2.  DISCUSSION 


The  most  recent  changes  to  PLOTEZ  have  been  made  during  the  period 
September,  1985  till  March,  1986.  The  author  had  started  working  on  the  MRL 
VAX-11/780  and  had  a  need  for  a  callable  plotting  routine  to  do  plots  of  ship 
magnetic  signatures.  By  using  a  set  of  magnetic  dipoles  to  represent  a 
ship's  magnetization  the  geomagnetic  field  perturbations  in  the  vicinity  of 
the  ship  could  be  calculated  [3]. 


The  code  PLOT  was  ported  to  the  MRL  VAX,  renamed  PLOTEZ  and  upgraded 
considerably  to  run  with  the  PLOTIO  IGL  library  [1],  The  PLOTIO  IGL  library 
provides  a  more  advanced  set  of  graphics  routines  than  the  earlier  PLOTIO  TCS 
library.  The  main  enhancements  implemented  in  this  new  version  are: 

(a)  software  generation  of  text  and  < b )  user  control  of  the  plot  aspect  ratio. 
Most  of  the  routines  have  been  thoroughly  rewritten  (in  FORTRAN  77)  so  that 
the  code  runs  faster,  is  more  efficient  and  better  documented  internally. 
Further  the  type  of  plot  output  (e.g  dashed  lines,  special  symbols  such  as 
stars)  has  been  greatly  expanded  using  the  capabilities  of  the  IGL  routines. 


3.  USING  PLOTEZ  ON  A  VAX  COMPUTER 


PLOTEZ  is  called  as  a  subroutine  from  a  FORTRAN  program  which 
contains  the  X,Y  values  stored  in  arrays.  For  example  the  following  test 
program  first  generates  sets  of  X,XE,Y,YE  values  ( XE, YE  are  the  errors  in  X,Y) 
and  then  calls  the  plotting  routine  for  output  to  a  graphics  terminal. 


PROGRAM  TEST 

INTEGER  NPT,  ICHAR,  IAX 

PARAMETER  (NPT=51) 

REAL  X(NPT) ,XE(NPT) ,Y(NPT) , YE(NPT) 

C 

C  GENERATE  THE  X,Y,XE,YE  VALUES 

DO  100  I  =  1 , NPT 
X(I)  =  FLOAT( 2*1-2) 

Y ( I )  =  100.*SIN(X(I)/3.0) 
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XE( I)  =  0.0 

YE{I)  =  0.1 *ABS(Y( I) ) 

1 00  CONTINUE 

C 

C  NOW  PLOT  THE  DATA 

TYPE  ICHAR, I AX  ?' 

ACCEPT  *, ICHAR, IAX 

CALL  PLOTEZ ( X, Y, XE, YE,NPT, ICHAR, IAX) 
C  .  STOP 

END 


In  the  call  to  PLOTEZ  the  arrays  X,Y,XE,YE  carry  the  data  and  NPT 
refers  to  the  number  of  points  of  the  arrays  to  be  plotted.  The  parameters 
ICHAR, IAX  specify  characteristics  of  the  plot  required  and  will  be  described 
in  the  next  section. 

For  a  user  running  on  the  MRL  VAX  the  program  TEST  must  now  be 
compiled  and  linked  with  PLOTEZ  and  the  IGL  library: 


$ FORTRAN  TEST 

SLINK  TEST, [RYANP. PLOT] PLOTEZ, SYS$LIBRARY : IGL/LIB 


To  run  the  program  it  is  necessary  to  be  logged  on  to  a  suitable 
graphics  terminal  and  type: 

$RUN  TEST 

For  a  user  at  a  different  site  it  would  be  necessary  to  first 
install  PLOTEZ  and  the  IGL  library,  then  compile  and  link  PLOTEZ  and  the  IGL 
library  with  the  user  program.  A  listing  of  the  program  is  given  in  an 
appendix  at  the  end  of  this  report. 


3.1  The  Parameters  ICHAR  and  IAX 


The  parameters  ICHAR  and  IAX  included  in  the  program  listing  above 
are  used  to  specify  the  type  of  plot  requested.  ICHAR  controls  the  plotting 
symbol  used.  There  are  3  main  ranges  of  integer  values  ICHAR  can  take: 

1.  0< ICHAR< 1 7  the  IGL  special  symbols  are  used.  For  example  for  ICHAR=4 
triangles  are  plotted  for  the  points;  for  ICHAR=1 2  stars  are  plotted 
for  the  points. 

2.  ICHAR=17,18  a  semi-  or  full  histogram  plot  is  done. 

3.  ICHAR>=20  the  points  are  joined  up  as  either  a  continuous  line 
(ICHAR=20)  or  a  dashed  line  (ICHAR>20). 

Note  that  for  ICHAR  <  0  error  bars  are  not  plotted. 
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The  parameter  I AX  specifies  the  graphic  action  requested.  Currently 
I AX  can  take  all  integer  values  from  -2  to  +4.  For  example,  if  IAX=-1  the 
program  asks  for  scaling  information,  titles  etc.  and  then  erases  the  screen 
drawing  a  new  frame  with  tick  marks  etc.  In  contrast,  for  IAX=1  the  screen 
is  not  erased,  the  same  scale  as  was  determined  by  the  last  plot  is  used  and 
the  data  over plotted  on  to  the  same  frame. 

Further  information  on  the  action  specified  by  the  parameters  ICHAR 
and  I AX  is  given  in  Appendix  A. 


3 . 2  Response  to  PLOTEZ  Prompts 

If  I AX  is  negative,  PLOTEZ  prompts  the  user  for  various  input 
parameters.  These  must  be  entered  (and  in  the  correct  format)  or  the  program 
will  wait  indefinitely  (or  crash  if  the  wrong  type  of  data  is  input). 


Question  1.  'PLOT  TYPET?  [P=AUTO  SCALE, S=SET  SCALE, E= EXIT, DEF=P] >> 


If  'P'  or  return  is  entered  the  program  scales  the  input  data 
automatical ly  whereas  if  'S'  is  entered  the  program  asks  for  the  minimum  and 
maximum  X,Y  values  to  be  plotted.  If  ' E'  is  entered  the  program  exits 
without  plotting.  The  single  character  input  here  can  be  in  either  upper  or 
lower  case. 


Question  2.  'FOR  SEMI-LOG  PLOT  TYPE  L  [ DEF=LINEAR] >> ' 

If  'L'  or  '1'  is  entered  the  Y  data  points  are  scaled 
logarithmically;  if  anything  else  is  entered  the  Y  data  points  are  scaled 
linearly. 


Question  3.  If  'S'  or  's'  has  been  entered  for  Q1  ,  then  the  program  asks:' 
MIN  X,  MAX  X,  MIN  Y,  MAX  Y  VALUES  TO  PLOT  [NO  DEF. ]  >> ' 

These  can  be  entered  in  free  format.  For  example  a  typical 
response  would  be  :  0.0,10.0,-100.0,100.0 


Questions  4,5,6. 


'TYPE  TOP  TITLE' >> 
'TYPE  X-AXIS  TITLE' >> 
'TYPE  Y-AXIS  TITLE' >> 


Here  the  3  titles  are  prompted  for  -  the  TOP  title  which  is  printed 
above  the  frame,  the  X-AXIS  title  printed  beneath  the  X-axis,  and  the  Y-AXIS 
title  which  is  printed  to  the  left  of  the  frame  and  rotated  through  90 
degrees.  After  each  title  is  entered  the  return  key  must  be  pressed. 
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The  code  will  accept  character  strings  of  up  to  80  characters  in 
length  as  input  to  these  prompts. 


Question  7.  'CHARACTER  SIZE  FOR  TOP,  X- ,  Y-AXIS  TITLES  [DEF=1 . , 1 . , 1 . ]  >>' 

The  size  of  the  characters  used  for  the  3  titles  can  be  specified 
separately  here.  If  return  is  entered  the  3  titles  are  printed  using  the 
default  character  size.  The  values  input  scale  the  title  sizes  from  the 
default.  For  example  if  2.0  is  entered  for  the  X-axis  character  size  and  1 .0 
for  the  Y-axis  character  size,  the  X-axis  title  characters  will  be  twice  as 
big  as  the  Y-axis  characters.  The  character  sizes  are  restricted  to  be  a 
maximum  value  of  2.0.  The  default  input  values  are  1.0, 1.0, 1.0  which  will  be 
set  if  the  return  key  is  pressed. 


Question  8.  'NUMBERS  OF  X-  AND  Y-AXIS  INCREMENTS  [DEF=5,5]' 

The  numbers  of  X-  and  Y-axis  increments  between  large  tick  marks  are 

entered  here.  If  NTX  or  NTY  exceeds  10  they  are  reset  to  10.  If  return  is 
entered  NTX  and  NTY  are  both  set  to  their  default  values  5.  A  typical 

response  to  this  prompt  would  be  :  6,10  specifying  6  increments  along  the  X- 

axis  and  10  along  the  Y-axis.  Note  that  if  log  scaling  is  requested  for  the 
Y-data  points  the-  value  of  NTY  will  be  calculated  automatically  and  the  value 
input  here  will  be  ignored. 


Question  9.  'X,Y  VIEWPORT  FACTORS  [DEF=1.,1.]' 

Here  the  factors  XVW.YVW  which  define  the  size  of  the  viewport  are 
entered.  The  viewport  is  the  area  on  the  display  (screen  or  plotter)  on  to 
which  the  plot  is  drawn.  A  typical  response  would  be  :  1.5, 1.2  specifying 
that  the  viewport  is  compressed  by  a  factor  of  1 . 5  in  the  X-direction  and  by 
1.2  in  the  Y-direction.  In  this  way  the  aspect  ratio  of  the  plot  can  be 
altered.  The  defaults  (entered  by  pressing  return)  are  both  equal  to  1 .0. 


Question  10.  'DEVICE,  OPTION  [ DEF=401 0, 1 ] ' 

Here  the  two  parameters  IDEV,I0PT  required  by  the  IGL  initialization 
routine  GRSTRT  (see  ref.  [1])  are  entered.  The  default  values  are  IDEV=4010, 
I0PT=1  (TEK401 0  or  compatible  terminal).  Note  that  if  hard  copy  output  is 
required  on  to  a  TEK4662  or  compatible  plotter  the  program  will  pause  after 
this  question  has  been  answered  to  give  the  user  time  to  set  up  the  plotter 
etc.  When  this  has  been  done  typing  any  integer  or  return  will  enable  the 
program  to  continue.  For  a  TEK4662  plotter  it  is  necessary  to  type  in  several 
more  returns  to  enable  it  to  start  plotting. 


4.  PROGRAM  DETAILS 


At  present  PLOTEZ  contains  about  900  lines  of  code  including  about 
350  lines  of  comments.  The  actual  code  contains  roughly  550  active  lines  of 
FORTRAN  77  and  its  incorporation  into  a  user  program  should  not  seriously  degrade 
performance.  The  maximum  number  of  points  which  can  be  plotted  in  a  single 
call  to  PLOTEZ  is  set  to  4096  using  a  PARAMETER  statement.  PLOTEZ  is  written 
to  conform  to  the  FORTRAN  77  standard  and  frequent  use  is  made  of  FORTRAN  77 
features  such  as  IF-THEN-ELSE-ENDIF  constructs,  CHARACTER  type  variables  etc. 


PLOTEZ  projects  the  plot  vectors  on  to  a  Tektronix  terminal  using 
the  default  system  of  GDUs  (Graphic  Display  Units)  as  the  display  surface 
units.  On  a  Tektronix  4010  terminal,  for  example,  the  default  window  in  GDUs 
is  roughly  130  units  along  the  X-axis  and  100  units  along  the  Y-axis.  The 
use  of  GDUs  means  that  the  output  is  independent  of  the  physical  size  of  the 
screen  since  each  GDU  is  1/1 00th  the  distance  of  the  shorter  axis  on  the 
display  surface. 


PLOTEZ  has  a  main  routine  PLOTEZ  and  6  subsidiary  routines.  The 
main  routine  accepts  the  input  data  from  the  calling  program,  decides  what 
type  of  plot  is  required  then  initializes  the  device  and  calls  the  subsidiary 
routines  to  do  the  scaling  and  plotting.  The  function  of  these  routines  is  as 
follows. 

(a)  FRAME  -  draws  the  frame  with  tick  marks  and  numbers  and  also  draws  the 
titles.  Note  that  since  the  text  is  software  generated  it  is  all  scaled  to 
the  same  size  as  the  plot. 

(b)  PLOTXY  -  plots  the  points  on  the  frame  drawn  by  FRAME.  The  type  of  plot 
symbol  is  specified  by  the  parameter  ICHAR. 

(c)  LOGTRAN  -  converts  the  Y-values  into  their  L0G1 0  representation  if  a 
semi-log  plot  is  requested. 

(d)  LINSCALE  -  determines  the  increments  between  the  tick  marks  and  the 
starting  position  on  each  axis. 

(e)  LOGSCALE  -  determines  increments  between  tick  marks  etc.  in  the  Y-axis 
for  semi-log  plots. 

(f)  XPCW  -  decomposes  the  input  real  variable  into  its  two  scientific- 
notation  components  and  returns  these  values  to  the  calling  routine.  For 
example,  if  -30.0  is  input  the  real  number  -3.0  and  the  integer  1  are 
returned.  This  is  used  by  the  scaling  routines. 

(g)  TICK  -  draws  tick  marks  or  line  segments  between  two  specified  points. 
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5.  SAMPLE  OUTPUT 


The  program  in  section  3  produces  the  output  shown  in  Fig.  1  when 
the  input  parameters  6,-2  are  specified  for  the  parameters  ICHAR,IAX.  The 
program  calculates  a  set  of  sinusoidal  Y-values  with  the  Y-error  bars 
proportional  to  the  absolute  Y-values.  The  program  plots  the  data  points  as 
diamonds  with  vertical  error  bars. 


A  more  complicated  example  is  given  in  Fig. 2  where  the  three-axis 
magnetic  field  perturbations  due  to  a  naval  vessel  passing  over  a  mine  located 
at  a  depth  of  20.0  m  are  plotted  on  the  same  graph.  The  X-field  component  is 
represented  by  the  continuous  line,  the  Y-field  component  by  the  dotted  line, 
and  the  Z-field  component  by  the  dash-dot  line.  The  ship  magnetization  is 
modelled  by  a  set  of  20  three-axis  magnetic  dipoles  spaced  6  m  apart  along  a 
midships  line  at  water  level.  The  X,Y,Z  field  components  are  calculated  by 
summing  the  contributions  from  each  dipole  in  each  axis  using  the  equations 
for  a  magnetic  dipole.  The  set  of  dipoles  used  is  a  test  set  of  data  given  in 
Table  3.1  of  [3]  to  examine  the  validity  of  the  technique  of  representing  a 
ship's  magnetization  by  a  set  of  discrete  three-axis  dipoles.  The  dipole 
field  equations  are  taken  from  the  program  listing  given  in  the  Appendix  of 
13]. 


This  plot  is  produced  by  3  calls  to  PLOTEZ.  The  first  call  uses 
the  ICHAR  value  of  20  to  produce  a  continuous  line  plot  of  the  X-field 
component  and  IAX=-1  to  set  the  scale  of  the  plot  manually.  Then  the  Y  and  Z 
components  of  the  field  are  overplotted  on  the  same  plot  by  calls  to  PLOTEZ 
with  ICHAR, IAX  values  of  21, +1  (dotted  line),  and  22, +1  (dash-dot  line).  The 
FORTRAN  code  fragment  which  produced  this  output  is: 

CALL  PLOTEZ ( XDISP , BX, ERRX , ERRY, NUMSTEP ,20,-1 ) 

CALL  PLOTEZ(XDISP, BY, ERRX, ERRY, NUMSTEP, 21 ,+1 ) 

CALL  PLOTEZ ( XDISP, BZ, ERRX, ERRY, NUMSTEP, 22, +1 ) 

The  arrays  XDISP, BX, BY, BZ, ERRX, ERRY  carry  the  X-values,  the  X,Y,Z 
field  components  and  the  errors  respectively,  while  NUMSTEP  specifies  the 
number  of  points  plotted  to  make  the  3  curves  on  the  plot. 


6.  CONCLUSION 


The  graphics  program  PLOTEZ  is  now  available  as  a  callable 
subroutine  on  the  MRL  VAX  11/780.  It  is  a  device-independent  program  and 
will  run  on  any  TEK-4010  compatible  terminal.  To  interface  PLOTEZ  to  a  user 
program,  it  is  necessary  to  link  the  program  with  PLOTEZ  and  the  IGL  PL0T1 0 
graphics  library  on  the  VAX  system. 
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APPENDIX  A 


A.  Action  Taken  By  Values  of  ICHAR  Parameter 

The  ICHAR  parameter  is  used  to  give  different  characters  for  the 
plotted  points: 

1.  For  ICHAR  <  17  one  of  the  characters  from  the  default  software  font 
with  Graphics  Text  Emulation  (3D)  (see  [1])  is  used.  These  are: 

ICHAR  CHARACTER 

1  Dot  or  bullet 

2  Square 

3  Octagon 

4  Triangle 

5  + 

6  Diamond 

7  Square  with  x-cross 

8  Lozenge 

9  Octagon  with  x-cross 

10  Square  with  x-cross 

1 1  Nabla  or  del 

1 2  Star 

13  Asterisk 

14  X-cross 

1  5  Op  arrow 

16  Down  arrow 


2.  For  ICHAR  =  17  a  semi -his tog ram  plot  is  requested;  for  ICHAR=18  a  full 
histogram  plot  is  drawn. 

3.  For  ICHAR  =  20  the  points  are  joined  as  a  continuous  line.  For 
ICHAR>20  one  of  the  possible  dashed  lines  is  drawn: 

ICHAR  LINE  TYPE 

21  . 

22  - 

23  - -  -  - 

24  _ _ _ _ 

25  -  -  -  - 

etc.  as  given  in  the  IGL  user  manual  [1]. 


B.  Action  Taken  By  IAX  Parameter 

I AX  specifies  which  type  of  plot  is  wanted. 

IAX  ACTION 

+4  No  plot  is  done;  only  the  scale  is  calculated 
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+3  Draws  the  frame  only  using  the  scale  input  from  the 

calling  program 

+2  Erases  the  screen,  then  using  the  same  scale  as  was 

calculated  previously  for  the  last  plot,  draws  the 
frame  and  plots  the  data 

+1  No  frame  drawn,  screen  not  erased;  data  overplotted  on  the 

frame  drawn'  by  last  plot 

0  Auto-scaled  plot 

-1  Prompts  for  scaling  parameters  etc;  recycles  after 

plotting 

-2  Same  as  -1  but  returns  after  plotting 
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APPENDIX  B 


SUBROUTINE  PLOTEZ ( X, Y, XE, YE, NPT, ICHAR, IAX) 

C 

C  PLOTTING  ROUTINE  FOR  TEKTRONIX-TYPE  TERMINALS  ON  VAX  11/780 

C  WHICH  USES  THE  PLOT1 0  IGL  LIBRARY. 

C 

C  X  —  X  VALUES 

C  Y  —  Y  VALUES 

C  XE  —  ERROR  IN  X  VALUES 

C  YE  —  ERROR  IN  Y  VALUES 

C  LB  —  BYTE  ARRAY  -  IF  TRUE  PLOT  POINT,  ELSE  NOT 

C  NPT  —  NUMBER  OF  POINTS 

C 


C 

ICHAR 

CONTROLS  PLOT  CHARACTER 

c 

0 

— 

ERROR  BARS  ONLY 

c 

1 

— 

POINTS 

c 

2 

— 

SQUARES 

c 

3 

— 

OCTAGONS 

c 

4 

— 

TRIANGLES 

c 

5 

— 

+  'S 

c 

6 

— 

DIAMONDS 

c 

7 

— 

SQUARE  WITH  X-CROSS 

c 

8 

— 

LOZENGE 

c 

9 

— 

OCTAGON  WITH  X-CROSS 

c 

10 

— 

SQUARE  WITH  4-CROSS 

c 

11 

— 

NABLA  OR  DEL 

c 

12 

— 

STAR 

c 

13 

— 

ASTERISK 

c 

14 

— 

X-CROSS 

c 

15 

— 

UP  ARROW 

c 

16 

— 

DOWN  ARROW 

c 

c 

r» 

— 

GIVES  SEMI-HISTOGRAM 

c 

18 

— 

GIVES  FULL  HISTOGRAM 

c 

c 

20 

— 

CONTINUOUS  LINE 

c 

>20 

— 

DASHED  LINES  : 

c 

21 

— 

DOTTED  LINE  { . ) 

c 

22 

— 

DASHED-DOT  LINE  (-•-•-•  ) 

c 

23 

— 

DASHED  LINE  ( - ) 

c 

24 

— 

LONG  DASHED  LINE 

c 

25 

— 

DASH  -  TWO  DOT  LINE  (-••- 

c 

ETC. 

SEE  IGL  MANUAL 

C 

C  IF  ICHAR  -VE  ERROR  BARS  ARE  NOT  PLOTTED 

C 


C  IAX  —  TYPE  OF  PLOT 

C  4  —  NO  PLOT,  CALCULATES  SCALE  ONLY 

C  3  —  PLOTS  AXES  ONLY.  SCALE  INPUT  FROM  CALL. 

C  2  —  PLOTS  AXES,  SAME  SCALE  AS  LAST  PLOT. 

C  1  —  NO  AXES,  NO  ERASE,  SAME  SCALE  AS  LAST  PLOT. 

C  0  — -  NORMAL  PLOT  (AUTO  SCALE, NO  TITLE  OPTIONS  ETC). 


B-1 


on  -*  n  n  o  o  o  no 


-1 


-2  — 


ASKS  FOR  PLOT,  WITH  OPTIONS  FOR  TITLE 
AND  USER  SCALING.  RECYCLES  AFTER  PLOTTING. 
ASKS  FOR  PLOT,  WITH  OPTIONS  FOR  TITLE 
AND  USER  SCALING.  RETURNS  AFTER  PLOTTING. 


C 
C 
C 
C 
C 
C 

C - 

C 
C 
C 
C 
C 
C 
C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

c  LATEST  VERSION  AS  OF  APR- 1986  P.RYAN 

C  PHYSICAL  CHEMISTRY  DIVISION 

C  MATERIALS  RESEARCH  LABORATORIES 

C  P.O.  BOX  50,  ASCOT  VALE  3032,  VICTORIA,  AUSTRALIA 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

C  SET  MAXIMUM  NUMBER  OF  POINTS  TO  PLOT 

PARAMETER  (MAXPTS  =  4096) 


FOR 

IAX  =  <3, 

4>  USE: 

XMIN 

—  MIN. 

X 

VALUE 

TO 

PLOT 

XMAX 

—  MAX. 

X 

VALUE 

TO 

PLOT 

YMIN 

—  MIN. 

Y 

VALUE 

TO 

PLOT 

YMAX 

—  MAX. 

Y 

VALUE 

TO 

PLOT 

DECLARATION  STATEMENTS 
DIMENSION  X(* ) ,  Y(*  ) ,XE(* ) ,  YE(* ) 

LOGICAL  LB( MAXPTS),  LOGPLOT,  SETSCALE 
REAL  LY( MAXPTS), LYE (MAXPTS) 

CHARACTER  *  80  TITLE, XTITLE, YTITLE 
CHARACTER  *  7  TOP,  X_AXIS,  Y_AXIS 
CHARACTER  *  1  ANS,  ANSLOG 

COMMON  /XYBLK/  XMIN, XMAX , YMIN, YMAX, XDELT, YDELT, XSIZE,  XWIDTH 
DATA  XLEFT , XRIGHT, YBOT , YTOP/ 4*0.0/ 

DATA  XVW,  YVW  /  1.0,  1.0/ 

DATA  IDEV,  IOPT/  4010,  1/ 

DATA  LOGPLOT,  SETSCALE  /.FALSE.,  .FALSE./ 

DATA  TOP,X_AXIS,Y_AXIS/'  TOP','  X-AXIS','  Y-AXIS'/ 

SET  DEFAULTS  : 

NTT, NXT,NYT  -  NO.  OF  CHARACTERS  IN  TOP , X- , Y-AXIS  TITLES 
NTX , NTY  -  NO.  OF  INCREMENTS  ALONG  X-,Y-AXES 

XSIZE, XWIDTH  -  CHARACTER  WIDTH  AND  CELL  WIDTH  ( =XSIZE+0. 385 ) 
00  NTT  =  0 
NXT  =  0 
NYT  =  0 
NTX  =  5 
NTY  =  5 
XSIZE  =  1 .410 
XWIDTH  =  1  .795 

IF( IAX  .EQ.  4)  SETSCALE  =  .FALSE. 

FDR  AUTO-SCALED  PLOT  OR  OVERPLOTTING  TRANSFER  CONTROL  TO  900 
IF(IAX.GE.O)  GO  TO  900 


no-*  no  on 


EXIT  IF  I AX  IS  OUT  OF  RANGE 
IF(IABS(IAX)  .GT.  4)  THEN 

TYPE  ICHAR  OUTSIDE  RANGE:  ICHAR  >  4  * 

GO  TO  10000 
END  IF 

PAUSE  IF  OUTPUT  DEVICE  IS  A  TEX4662  PLOTTER 
IF( IDEV.EQ.4662)  ACCEPT  150,NUM 
50  FORMAT ( II  0) 

IAX  IS  NBG....  ASK  FOR  GRAFIX  PARAMS 
TYPE  200 

200  FORMAT (  '  PLOT  TYPE  ? ' , 

+  '  [P=AUTO  SCALE,  S=SET  SCALE,  E=EXIT,  DEF=P]  >>',$) 

ACCEPT  800, ANS 

IF(ANS.EQ. 'E* .OR.ANS.EQ. '  N' .OR.ANS.EQ. ' e ' .OR.ANS.EQ. 'n' ) 
+  GO  TO  10000 

IF(ANS.NE. ' E' .OR.ANS.NE. 'N' .OR.ANS.EQ. ’S' .OR. 

+  ANS.ne. 'e' .or.ANS.ne. 'n' .or.ANS.eq. ' s' )GO  TO  250 
GO  TO  10000 
250  TYPE  300 

300  FORMAT ( '  FOR  SEMI-LOG  PLOT  TYPE  L  [DEF=L INEAR]  >>',$) 

ACCEPT  800, ANSLOG 
LOGPLOT  =  .FALSE. 

IF(ANSLOG.EQ. 'L' .OR. ANSLOG. EQ. ' 1* )  LOGPLOT  =  .TRUE. 

IF  (ANS. EQ. 'S'. OR.ANS.EQ. 'S')  THEN 
SETSCALE  *  .TRUE. 

TYPE  350 

350  FORMAT ( $ '  MIN  X,  MAX  X,  MIN  Y,  MAX  Y  VALUES  TO  PLOT', 

+  '  [NO  DEF.]  »  ') 

ACCEPT  *,XMIN,XMAX, YMIN, YMAX 
ELSE 

SETSCALE  =  .FALSE. 

END  IF 

C 

C  INPUT  TITLING  INFORMATION 

C  NTT, NXT, NYT  -  NO  OF  CHARACTERS  IN  TITLES 

C  TITLE,  XTITLE,  YTITLE  -  TITLE  CHARACTER  STRINGS 

C  CHTOP,CHAXIS,CHYAXIS  -  TITLE  CHARACTER  SIZES 

TYPE  4 00, TOP 

400  FORMAT ( '  TYPE  ',A,'  TITLE  >>',$) 

ACCEPT  850, NTT, TITLE 
TYPE  400, X_AXIS 
ACCEPT  850, NXT, XTITLE 
TYPE  400, Y_AXIS 
ACCEPT  850, NYT, YTITLE 
TYPE  450 

450  FORMAT ( '  CHARACTER  SIZES  FOR  TOP,  X-, Y-AXIS  TITLES  ', 

+  '  [DEF=1 • , 1 . , 1 . ]  >>’,$) 

ACCEPT  500,CHTOP,CHXAXIS,CHYAXIS 
500  FORMAT( 3G20.7) 

CHTOP  =  MIN ( 2.0, CHTOP ) 

CHXAXIS  =  MIN(2.0,CHXAXIS) 

CHYAXIS  =  MIN(2.0,CHYAXIS) 
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550 

600 


650 

700 


750 


C 

c 

800 

850 

900 


950 


C 

C 

c 

1  000 


IF(CHTOP.LE.O.O)  CHTOP  =  1.0 
IF(CHXAXIS.LE.O.O)  CHXAXIS  =  1.0 
IF(CHYAXIS.LE.O.O)  CHYAXIS  =  1.0 
TYPE  550 

FORMAT(  $ , 1  NUMBERS  OF  X  AND  Y  AXIS  INCREMENTS  [DEF=  5,5]  >>',$) 

ACCEPT  600,NTX,NTY 

FORMAT(2I10) 

NTX  =  MIN( 1 0,NTX) 

NTY  =  MIN( 1 5, NTY) 

IF(NTX.LE.O)  NTX  =  5 
IF(NTY.LE.O)  NTY  =  5 
TYPE  650 

FORMAT ( '  X, Y  VIEWPORT  FACTORS  [DEF=1.,1.]  >>  ',$) 

ACCEPT  700,XVW, YVW 
EORMAT( 2G20. 7 ) 

XVW  *  MAX ( 1 . 0 , XVW ) 

YVW  =  MAX ( 1 .0, YVW) 

TYPE  750 

EORMAT( '  DEVICE,  OPTION  [DEF=4010,1]  >>',$) 

ACCEPT  600, IDEV, IOPT 
IF(IDEV.LE.O)  IDEV  =  4010 
IF( IOPT.LE. 0 )  IOPT  =  1 


FORMATS  USED  FOR  CHARACTER  DATA  INPUT 
FORMAT (A1) 

FORMAT (  Q,  A) 

IF  (SETSCALE  .E Q.  .TRUE.)  GO  TO  1100 
IF  (SETSCALE  .EQ.  .FALSE.)  THEN 
IF( IAX.GE. 1  .AND.  IAX.LE.3)  THEN 
TYPE  950 

FO  RMAT ( 1  *  *  *  *  WARNING  -  SCALES  NOT  YET  SET  SO  RETURN  ****') 
GO  TO  10000 
ELSE 

GO  TO  (1150,1150,1100,1000)  IAX 
END  IF 
END  IF 


SET  SCALE  BY  DETERMINING  MAX,  MIN  X,Y  VALUES  TO  PLOT 

OONTINUE 

XMIN  =  X(  1  ) 

XMAX  =  X(  1  ) 

YMIN  =  Y( 1 ) 

YMAX  =  Y(  1  ) 

DO  1050  1=1, NPT 
IF  ( ICHAR  .LT.  0)  THEN 
XMIN  =  MIN ( XMIN, X( I ) ) 

XMAX  =  MAX ( XMAX , X ( I ) ) 

YMIN  =  MIN ( YMIN, Y( I ) ) 

YMAX  =  MAX ( YMAX , Y ( I ) ) 

ELSE  IF  (ICHAR  .GE.  0)  THEN 
XMIN  =  MIN( XMIN, X(I)-XE(I)) 

XMAX  =  MAX ( XMAX , X ( I ) +XE ( I ) ) 

YMIN  =  MIN(YMIN,Y(I) -ABS ( YE( I ) ) ) 

YMAX  =  MAX ( YMAX , Y ( I ) + ABS ( YE ( I ) ) ) 
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END  IF 

1 050  CONTINUE 

SETSCALE  =  .TRUE. 

C 

C  LINEAR  SCALE  FOR  X-VALUES 

1100  CALL  LINSCALE(XMIN,XMAX, XMV , XDIV, IPX,NTX ) 

C 

C  LINEAR  OR  LOG  SCALE  FOR  Y-VALUES 

IF (LOG PLOT  .EQ.  .FALSE.)  CALL  LINSCALE( YMIN, YMAX, YMV, YDIV, IPY,NTY) 
IF( LOGPLOT  .EQ.  .TRUE.  )  CALL 
LOGSCALE ( YMIN , YM AX , YMV, NTY , IP Y , INC , ISPEC ) 

C 

C  TOR  LOG  SCALING  CONVERT  Y-VALUES  TO  LOG1 0  BASE 

1150  IF (LOGPLOT  .EQ.  .TRUE.  )  CALL  LOGTRAN (X,Y,XE,YE,LY,LYE, NPT ) 

C 

C  RETURN  IF  IAX  =  4 

IF( I AX. EQ. 4) GO  TO  10000 
C 

C  INITIALIZE  DEVICE  -  SET  REQUIRED  DEVICE  DRIVER 

C  ALSO  ENSURE  SOFTWARE  PRODUCED  TEXT  AND  SET  TEXT  SIZE 

C 

1200  IF ( IDEV.EQ.4662)  ACCEPT  1250,NUM 

1250  FORMAT (110) 

CALL  GRSTRT( IDEV, IOPT) 

CALL  TXQUAL ( 4 ) 

IF  ( IAX  .NE.  1  )  THEN 
XSIZE  =  XSIZE/XVW 
END  IF 

CALL  TXSIZE(0,  XSIZE,  0.0) 

CALL  BELL 

IF( IAX.NE. 1 )  CALL  NEWPAG 
IF(IAX.EQ.I)  GO  TO  1350 
C 

C  SET  WINDOW  FOR  PLOTTING 

C  XL EFT, XRIGHT, YBOT, YTOP  CORRESPOND  TO  THE  LIMITS  OF  THE  X,  Y 

C  VALUES  WHICH  CAN  BE  REFERENCED  ON  THE  SCREEN.  THESE  ARE  SET 

C  TO  ALLOW  AMPLE  ROOM  TO  WRITE  TITLES,  NUMBERS  ON  THE  AXES  ETC. 

C 

1300  XDELT  =  ABS ( XMAX-XMIN) 

XLEFT  =  XMIN  -  0.20*XDELT 
XRIGHT=  XMAX  +  0.10*XDELT 
IF (.NOT.  LOGPLOT)  THEN 
YDELT  *  ABS (YMAX- YMIN) 

YBOT  =  YMIN  -  0.1 5* YDELT 
YTOP  =  YMAX  +  0. 1 5* YDELT 
ELSE  IF (LOGPLOT)  THEN 

YDELT  =  ABS ( LOG1 0( YMAX ) -LOG1 0( YMIN) ) 

YBOT  =  LOGIO(YMIN)  -  0.15*YDELT 
YTOP  =  LOGIO(YMAX)  +  0.1 5* YDELT 
YMIN  =  LOGIO(YMIN) 

YMAX  =  LOG1 0 ( YMAX ) 

END  IF 

1350  CALL  WINDOW(XLEFT, XRIGHT, YBOT, YTOP) 

C 

C  XGDUMAX,  YGDUMAX  CORRESPOND  TO  THE  LIMITS  OF  THE  DISPLAY  SURFACE  IN 


C  GDUS  WHERE  THE  PLOT  IS  TO  BE  DRAWN.  THE  DISPLAY  SURFACE  IS  ROUGHLY 

C  130  GDUS  IN  THE  X-AXIS  AND  100  IN  THE  Y-AXIS. 

XGDUMAX  =  130./XVW 
YGDUMAX  =  1 OO./YVW 

CALL  VWPORT( 0.0, XGDUMAX ,0.0, YGDUMAX) 

C 

C  CALL  SUBROUTINE  PLOTXY  TO  PLOT  POINTS 

C 

1400  IF ( IAX.NE. 3  )  THEN 

IF  (.NOT.  LOGPLOT)  CALL  PLOTXY( X, Y,  XE, YE, LB ,NPT, ICHAR, LOGPLOT) 

IF  (LOGPLOT)  CALL  PLOTXY(X, LY,XE, LYE, LB, NPT, ICHAR, LOGPLOT) 

END  IF 
C 

C  CALL  SUBROUTINE  FRAME  TO  DRAW  AXES  AND  NOS.  ON  AXES. 

C 

1450  IF( IAX.NE.  1  ) CALL  FRAME ( XMV, YMV, XDIV, YD IV, IPX, IP Y, NTX, NTY, 

+  LOGPLOT, INC, ISPEC, 

+  TITLE, XTITLE, YTITLE,NTT, NXT, NYT,CHTOP ,CHXAXIS,CHYAXIS) 

C 

C  RETURN  ONCE  FRAME  ETC.  DRAWN  FOR  IAX=3 

IF( IAX  .EQ.  3)  GO  TO  9000 
C 

1500  IF( IAX.EQ.-1 )THEN 

CALL  HOME 
CALL  GRSTOP 
GO  TO  100 
END  IF 
C 

9000  CONTINUE 

CALL  HOME 
CALL  GRSTOP 
C 

10000  CONTINUE 

RETURN 
END 

Q********************************************************************** 

c 

C  SUBROUTINE  FRAME 

C 

C  DRAWS  AXES,  TICK  MARKS  AND  TEXT 

C* ********************************************************************* 

SUBROUTINE  FRAME ( XMV, YMV, XDIV, YDIV, IPX, IP Y, NTX, NTY, 

+  LOGPLOT, INC, ISPEC, TITLE, XTITLE, YTITLE,  NTT, NXT,NYT, 

+  CHTOP,CHXAXIS,CHYAXIS) 

C 

DIMENSION  YLOG( 20) ,  YMARK( 20 ) 

REAL  XVALUE( 20 ) , YVALUE( 20) ,XPOSX( 20) , YPOSY ( 20 ) 

INTEGER  NBEFOREX ( 20 ) ,  NBEFOREY ( 20 ) 

CHARACTER  *  80  TITLE,  XTITLE,  YTITLE 
LOGICAL  LOGPLOT 
REAL  XDEL,  YDEL 

COMMON  /XYBLK/  XMIN, XMAX , YMIN, YMAX , XDELT, YDELT, XS IZE,  XWIDTH 
C 

C  FIRST  DRAW  FRAME 


o  o 


CALL  SKIP 

CALL  MOVE( XMIN, YMIN) 

CALL  DRAW(XMAX, YMIN) 

CALL  DRAW ( XMAX , YMAX ) 

CALL  DRAW (XMIN, YMAX) 

CALL  DRAW (XMIN, YMIN) 

CALCULATE  X,  Y  INCREMENT  UNITS  CONVENIENTLY  SCALED. 

XDEL  =  XDELT/1 00. 

YDEL  =  YDELT/1 00. 

C 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
C  DRAW  X-AXIS  TICK  MARKS  +  NOS 

C  NOTE  :  XV  IS  THE  ACTUAL  X-VALUE 

C  XVALUE( )  IS  THE  SCALED  X-VALUE 

C  XINC  IS  THE  X- INCREMENT  BETWEEN  LARGE  TICK  MARKS 

C  XDIV  IS  THE  SCALED  INCREMENT  BETWEEN  LARGE  TICK  MARKS 

C  XMV  IS  THE  SCALED  MINIMUM  VALUE 

C  IPX,  IPXM  ARE  EXPONENTS  WHICH  SCALE  THE  X-VALUES 

C 

IF( IPX.LT.O)  IPXM  =  IPX 
IF( IPX.GE.O. AND. IPX.LE. 3 )  IPXM=0 
IF( IPX.GE.4)  IPXM  =  IPX 
XINC  =  XDIV  *  10.**IPX 

C 

C  FIRST  DRAW  X-AXIS  TICK  MARKS 

DO  100  1=1 , NTX+1 

XV  = ( XMV+XD IV  *  FLOAT (1-1  ))*10.**IPX 
IF  (I.GT.1  .AND.  I.LT.NTX+1)  THEN 
CALL  TICK (XV, YMIN, XV, YMIN+3.0* YDEL) 

CALL  TICK ( XV , YMAX- 3 . 0* YDEL, XV , YMAX ) 

END  IF 

IF  (I.LE.NTX)  THEN 

CALL  TICK( XV+0. 5*XINC, YMIN, XV+0. 5*XINC , YMIN+1 . 5* YDEL ) 

CALL  TICK (XV+0. 5* XINC, YMAX-1  . 5* YDEL, XV+0 . 5* XINC, YMAX) 

END  IF 

1 00  CONTINUE 

C 

C  NOW  SCAN  X-VALUES  TO  DETERMINE  THE  MAX.  NO  OF  SIGNIFICANT  FIGS 

C  REQUIRED  (NSIGX)  AND  STORE  THE  XVALUES,  THE  STARTING  POSITIONS  AND 

C  THE  NO.  OF  FIGS  BEFORE  THE  DECIMAL  PT.  IN  ARRAYS  XVALUE( ) ,  XPOSX(), 

C  AND  NBEFOREX  RESPECTIVELY.  NOTE  X-VALUES  ROUNDED  TO  2  DECIMAL  PLACES 

NSIGX  =  -1 
DO  1 25  I  =  1 , NTX+1 

XV  =  (XMV+XD IV* FLOAT ( 1-1  ) ) *1  0.0** IPX 
XVAL  =  XV  /  10.** IPXM 

XVAL  =  FLOAT (NINT(  1  00.0*XVAL) )  /  100.0 
C 

C  DETERMINE  HOW  MANY  FIGS.  AFTER  DECIMAL  PT. 

C  NFRAC  =  FRACTIONAL  PART  OF  XVAL  *  100 

C  IDIGIT1 ,  IDIGIT2  =  DIGITS  OF  NFRAC 

C  NBEFOREX  =  NO.  OF  FIGS  BEFORE  DECIMAL  PT. 

C  NAFTERX  =  NO.  OF  FIGS  AFTER  DECIMAL  PT. 

NFRAC  =  ABS(NINT( 1 00.* XVAL) -10 0*1 NT (XVAL) ) 

IDIGIT1  =  NFRAC  /  10 


on-*  on-*  oooo 


IDIGIT2  =  NFRAC  -  10+IDIGIT1 
NAFTERX  =  -1 

IF  (IDIGIT1  .GT.  0)  NAFTERX  =  1 
IF  ( IDIGIT2  .GT.  0)  NAFTERX  =  2 
NSIGX  =  MAX (NAFTERX, NSIGX) 

DETERMINE  START  OF  CHARACTER  STRING  USING  KNOWLEDGE  THAT  CHARACTER 

ELL 

IS  XWIDTH*XDEL  GDU'S  WIDE,  THEN  STORE  NBEPOREX( ) , XPOSX( ) , XVALUE( ) 
CALL  XPOW ( XVAL , BX , IPXVAL ) 

NBEFOREX ( I )  =  MAX ( 1 , 1 +IPXVAL )  +  (1-SIGN( 1 .0,XVAL) )/2 
XPOSX(I)  =  XV  -  FLOAT ( NBEFOREX ( I ) ) *XWIDTH*XDEL 
XVALUE( I)  =  XVAL 
25  CONTINUE 

NOW  PRINT  X- VALUES 
DO  150  I  =  1 , NTX+1 
CALL  MOVE(XPOSX(I) ,YMIN-4.0*YDEL) 

NCHARX  =  NBEFOREX ( I)  +  NSIGX  +  1 
CALL  RNUMBR( XVALUE( I ) , NSIGX, NCHARX ) 

50  CONTINUE 

WRITE  10  RAISED  TO  POWER  IPXM  IF  IPXM  IS  NON-ZERO 
IF  ( IPXM.NE. 0 )  THEN 

CALL  MOVE  (XMAX-2.*XDEL,YMIN-1 2.0*YDEL) 

CALL  TEXT( 3, ' XI  O' ) 

CALL  MOVE(XMAX+1 . 0*XDEL, YMIN-9. 0* YDEL ) 

CALL  INUMBR( IPXM, 3) 

END  IF 

C 

CYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYY 
C  NOW  DRAW  Y-AXIS  TICK  MARKS  AND  NOS. 

C  SEE  COMMENTS  ABOVE  FOR  X-TICK  MARKS  ETC.  THAT  IS 

C  LOOP  200  DRAWS  TICK  MARKS,  LOOP  225  SCANS  THE  Y-VALUES  TO 

C  DETERMINE  THE  NO.  OF  SIG.  FIGS  REQUIRED,  AND  THEN  LOOP  250 

C  PRINTS  THE  Y-VALUES 

C 

IF( LOGPLOT)  GO  TO  300 
IF(IPY.LT.O)  IPYM  =  IPY 
IF( IP Y.GE.O. AND. IPY.LE.4)  IPYM  =  0 
IF(IPY.GE.S)  IPYM  =  IPY 
YINC  =  YDIV  *  1 0. **IPY 
C 

C  FIRST  DRAW  TICK  MARKS 

DO  200  J=1 , NTY+1 

YV  =  ( YMV+YDIV*FLOAT( J-1 ) )*10.**IPY 
IF  (J.GT.1  .AND.  J. LT. NTY+1 )  THEN 
CALL  TICK( XMIN, YV, XMIN+2 • *XDEL, YV ) 

CALL  TICK( XMAX-2 . *XDEL, YV, XMAX , YV) 

END  IF 

IF  ( J.LE.NTY)  THEN 

CALL  TICK( XMIN, YV+0. 5* YINC, XMIN+XDEL,  YV+0.5*YINC) 

CALL  TICK(XMAX-XDEL, YV+0, 5* YINC, XMAX, YV+0. 5* YINC) 

END  IF 

200  CONTINUE 
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C  NOW  SCAN  Y- VALUES  ROUNDED  TO  2  DECIMAL  PLACES.  SEE  COMMENTS  IN  LOOP 

C  125  ABOVE 

NSIGY  =  -1 
DO  225  J  =  1 , NTY+1 

YV  =  (YMV+YDIV*FLOAT( J-1 ) )*10.0**IPY 
YVAL  =  YV  /  1 0. **IPYM 

YVAL  =  FLOAT(NINT( YVAL*1 00.0) )  /  100.0 
NFRAC  =  ABS( NINT( 1 00. *YVAL)-1 00*INT( YVAL) ) 

IDIGIT1  =  NFRAC  /  10 
IDIGIT2=  NFRAC  -  10*IDIGIT1 
NAFTERY  *  -1 

IF  ( IDIGIT1  .GT.  0)  NAFTERY  =  1 
IF  ( IDIGIT2  .GT.  0)  NAFTERY  =  2 
NSIGY  =  MAX (NAFTERY, NSIGY) 

CALL  XPOW ( YVAL , B Y , IP YVAL ) 

NBEPOREY(J)  =  MAX( 1 , IPYVAL+1 )  +  ( 1 -SIGN( 1 . 0 , YVAL ) ) /2 
YPOSY(J)  =  YV 
YVALUE(J)  =  YVAL 
225  CONTINUE 

C 

C  NOW  PRINT  Y- VALUES 

DO  250  J  =  1, NTY+1 
NOTARY  =  NBEFOREY ( J )  +  1  +  NSIGY 
XPOS  =  XMIN  -  (FLOAT(NCHARY)+0.5)*XWIDTH*XDEL 
CALL  MOVE ( XPOS , YPOSY ( J ) ) 

CALL  RNUMBR( YVALUE( J) , NSIGY, NCHARY) 

250  CONTINUE 

C 

C  WRITE  10  RAISED  TO  THE  POWER  IPYM  IF  IPYM  NON-ZERO 

IF  ( IPYM. NE. 0 )  THEN 

CALL  MOVE( XMIN-1 0. *XDEL, YMAX+4. 0*YDEL ) 

CALL  TEXT( 3 , 'X10' ) 

CALL  MOVE( XMIN-7 .0*XDEL, YMAX+6. 0* YDEL ) 

CALL  INUMBR( IPYM, 3 ) 

END  IF 
GO  TO  9999 
C 

CLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLL 
C  DRAW  Y-AXIS  TICK  MARKS  AND  NOS  FOR  LOG  SCALING 

C 

300  CONTINUE 
C 

IF( ISPEC.EQ.O)  THEN 

C  DRAW  TICK  MARKS,  NOS  FOR  ISPEC  =  0  I.E.  IN  POWERS  OF  10 

DO  400  I LOG  =  1, NTY+1 
IVAL  =  IPY  +  ( ILOG-1 ) *INC 
YVALOG  =  FLO AT (IVAL) 

CALL  MOVE( XMIN-8. *XDEL, YVALOG) 

CALL  TEXT( 2 , '10' ) 

CALL  MOVE( XMIN-6. 5*XDEL, YVALOG+3. 5* YDEL ) 

CALL  INUMBR(IVAL, 3) 

IF(IL0G.GT.1  .AND.  ILOG. LT. NTY+1 )  THEN 

CALL  TICK (XMIN, YVALOG, XMIN+2 . *XDEL , YVALOG) 

CALL  TICK ( XMAX , YVALOG , XMAX- 2 . *XDEL, YVALOG) 
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END  IF 

IF( ILOG.LT. NTY+1  )  THEN 
MARK  =  (2*INC-3)/2  +  1 
IF( INC . EQ. 1 )  YFAC  =5.0 
IF( INC.NE. 1 )  YFAC  =  10.0 
DO  350  IMARK  =  1 , MARK 
YVALOG  =  IVAL  +  LOCI 0( ( YFAC** IMARK) ) 

CALL  TICK( XMIN, YVALOG, XMIN+1 .0*XDEL, YVALOG) 

CALL  TICK(XMAX, YVALOG, XMAX-1 . 0* XDEL , YVALOG ) 

CONTINUE 
END  IF 
CONTINUE 

FOR  CASE  OF  ISPEC  =  1 

FIRST  CALCULATE  YLOG ( ) ,  YMARK( )  VALUES  WHERE 
YLOG  -  THE  LOGS  OF  THE  Y-LOCATIONS  OF  THE  MAIN  TICK  MARKS 
YMARK  -  THE  VALUE  1.,2.,,  OR  5  WRITTEN  NEXT  TO  THE  TICK  MARK 
ELSE  IF  ( ISPEC  .EQ.  1  )  THEN 
YLOG( 1 )  =  YMIN 
YLOG( NTY+1 )  =  YMAX 
YMARK ( 1 )  =  YMV 
DO  450  IMARK=2, NTY+1 
YMARK (IMARK)  = YMARK ( IMARK- 1 ) *  2 . 0 

I F ( YMARK ( IMARK ). NE . 2 .. AND . YMARK ( IMARK ). NE . 1 0 .) YMARK ( IMARK )= 5 . 

IF (YMARK (IMARK)  .EQ.  10.)THEN 
YMARK ( IMARK )  =  1.0 
IPY  =IPY  +  1 
END  IF 

YTEMP  *  YMARK( IMARK)  *  10.**IPY 
YLOG (IMARK)  =  LOG10( YTEMP) 

CONTINUE 

NOW  DRAW  TICK  MARKS 
DO  550  JMARK  =  1 , NTY 

IF( YMARK (JMARK) . EQ. 1 .0  .AND.  JMARK. GT. 1  .AND.  JMARK. LT. NTY+1 )  THEN 
CALL  TICK ( XMIN, YLOG ( JMARK) ,XMIN+2.0*XDEL, YLOG (JMARK) ) 

CALL  TICK ( XMAX, YLOG ( JMARK ), XMAX- 2. 0*XD EL, YLOG (JMARK) ) 

ELSE  IF( YMARK (JMARK)  .NE.  1.0  .AND.  YLOG ( JMARK )  .LT.  YMAX)  THEN 
MARK1  =  NINT( YMARK (JMARK ) ) 

MARK 2  =  2* MARK 1  -  MOD ( MARK 1,2) 

DO  500  IMAR  =  MARK1  , MARK2 

YVALOG  =  LOG1 0( FLOAT ( IMAR )/ FL0AT(MA RK1 ) )  +  YLOG( JMARK) 

CALL  TICK (XMIN, YVALOG, XMIN+XDEL, YVALOG) 

CALL  TICK( XMAX , YVALOG , XMAX-XDEL , YVALOG ) 

CONTINUE 
END  IF 
CONTINUE 

NOW  DRAW  NUMBERS  ON  Y-AXIS  FOR  LOG  SCALING  FOR  ISPEC  =1 
DO  600  KMARK=1 , NTY+1 

CALL  XPOW (10. **YLOG( KMARK) , YMANT, I  EXP) 

MANTISS  =  NI NT (YMANT) 

IF  ( MANTISS. NE. 1 )  THEN 
IF(NTY.LE.IO)  THEN 

CALL  MOVE( -3. 0*XWIDTH*XDEL+XMIN, YLOG ( KMARK ) ) 


0 


CALL  INUMBR(MANTISS, 2) 

END  IF 

ELSE  IF  (MANTISS.EQ. 1 )  THEN 

CALL  MOVE( XMIN-3 • 5*XWIDTH*XDEL, YLOG (KMARK) ) 

CALL  TEXT( 2, '10') 

CALL  MOVE( XMIN-2. 5*XWIDTH*XDEL, YLOG ( KMARK )+ 3 . 0* YDEL ) 

CALL  INUMBR(IEXP,2) 

END  IF 

600  CONTINUE 

C 

C  SPECIAL  CASE  FOR  SCALING  FROM  2.*10.**IPY  TO  5.*10.**IPY 

IFCNTY.EQ.1  .AND.  YMV.EQ.2.0)  THEN 
IF ( IPY  .NE.  0)  THEN 

CALL  MOVE( XMIN-1 4. *XDEL, YMAX-1 0.*YDEL) 

CALL  TEXT( 3, 'XI  O' ) 

CALL  MOVE ( XMIN-1 0. *XDEL , YMAX-7 . *  YDEL ) 

CALL  INUMBR( IPY, 2) 

END  IF 

DO  650  IP  =  3,4 

YIP  =  LOG1 0( FLOAT ( IP) )  +  IPY 

CALL  MOVE( XMIN-3. 5*XWIDTH*XDEL, YIP) 

CALL  INUMBR( IP , 2 ) 

650  CONTINUE 

END  IF 
C 

END  IF 

9999  CONTINUE 
C 

CTTTTTTTTTTOTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTT 
C  DRAW  TITLES 

C 

C  MAIN  TITLE 

IF(NTT.GT.O)  THEN 
CALL  TXANGUO.  ) 

CALL  TXSIZE( 0 , CHTOP*XSIZE, 0.0) 

CALL  MOVE( XMIN+5 . 0*XDEL, YMAX+2 . 0*YDEL ) 

CALL  TEXT( NTT, TITLE) 

END  IF 

C 

C  X-TITLE  -  TITLE  BELOW  X-AXIS 

IF(NXT.GT.O)  THEN 

CALL  MOVE( XMIN+20. *XDEL, YMIN-1 2.5*YDEL) 

CALL  TXANGL ( 0 . ) 

CALL  TXS IZE( 0, XSIZE*CHXAXIS, 0.0) 

CALL  TEXT( NXT, XTITLE) 

END  IF 
C 

C  Y-TITLE  -  TITLE  TO  LEFT  OF  Y-AXIS  ROTATED  THROUGH  90  DEG. 

IF(NYT.GT.O)  THEN 
CALL  TXANGL (90.) 

IF( LOGPLOT)  CALL  MOVE( XMIN-1 7 . 5*XDEL, YMIN+1 5 . *YDEL ) 

IF( .NOT.  LOGPLOT)  THEN 

XLOC  =  XMIN  -  FLOAT (  IP Y-IP YM+NSIGY+6 ) *XDEIL* XWIDTH 
CALL  MOVE( XLOC, YMIN+1 5.* YDEL) 

END  IF 
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CALL  TXS IZE{ 0 , XSIZE*CHYAXIS . 0.0) 

CALL  TEXT( NYT, YTITLE) 

END  IF 
C 

RETURN 

END 

Q********************************************************************** 

C  SUBROUTINE  PLOTXY 

C 

C  PLOTS  THE  (  X ,  Y)  POINTS  ON  THE  FRAME  DRAWN  BY  FRAME 

C 

^********************************************************************** 

SUBROUTINE  PLOTXY( X, Y, XE, YE, LB, NPT, ICHAR , LOGPLOT) 

REAL  X( * ) , Y( * ) , XE( * ) , YE( * ) 

LOGICAL  LB( * ) ,  LOGPLOT 

COMMON  /XYBLK/  XMIN, XMAX, YMIN, YMAX, XDELT,  YDELT,  XSIZE, XWIDTH 

TEST  WHICH  POINTS  LIE  WITHIN  RANGE  AND  SET  LB  VALUES  TO 
TRUE  OR  FALSE  FOR  EACH  (X,Y)  PAIR 

DO  50  IP  =  1  ,  NPT 
LB(IP)  =  .FALSE. 

3  CONTINUE 

DO  1 00  JP  =  1 , NPT 

FOR  ICHAR  <  0  EXCLUDE  ERROR  BARS  FROM  CONSIDERATION 
IF  (ICHAR  .LT.  0)  THEN 

IF  ( X ( JP ) . GE. XMIN  .AND.  X ( JP) .LE.XMAX  .AND. 

1  Y ( JP ) .GE. YMIN  .AND.  Y ( JP ) . LE. YMAX ) 

2  LB(JP)  =  .TRUE. 

FOR  ICHAR  GE  0  AND  LINEAR  Y- SCALE  INCLUDE  SIZE  OF  X,Y  ERROR  BARS 
ELSE  IF  (ICHAR  .GE.  0  .AND.  LOGPLOT  .EQ.  .FALSE.)  THEN 

IF  ( X( JP)-XE( JP) .GE.XMIN  .AND.  X ( JP ) +XE ( JP ) . LE. XMAX  .AND. 

1  Y( JP)-YE( JP) .GE.YMIN  .AND.  Y ( JP )+ YE ( JP ) . LE. YMAX ) 

2  LB(JP)  =  .TRUE. 

FOR  ICHAR  GE  0  AND  LOG  Y-SCALE  INCLUDE  SIZE  OF  X-ERROR  BARS  ONLY 
ELSE  IF  (ICHAR  .GE.  0  .AND.  LOGPLOT  . EQ.  .TRUE.)  THEN 

IF  (  X(JP)-XE(JP)  .GE.XMIN  .AND.  X ( JP ) +XE ( JP ) . LE. XMAX  .AND. 

1  Y(JP) .GE.YMIN  .AND.  Y(JP)  .LE. YMAX ) 

2  LB(JP)  =  .TRUE. 

END  IF 

00  CONTINUE 

TEST  WHICH  TYPE  OF  PLOT  IS  REQUIRED 
ITYPE  =  ABS( ICHAR) 

I F( I TYPE. EQ. 0 )  GO  TO  550 

IF( ITYPE. GT. 0  .AND.  ITYPE. LE. 16)  GO  TO  150 
IF( ITYPE. GE.  1  7. AND.  ITYPE. LE. 18)  GO  TO  250 
IF( ITYPE. GE. 20  )  GO  TO  350 

C  PLOT  SPECIAL  SYMBOLS  USING  THE  MARKER  IGL  ROUTINE 
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C  TOR  THE  SYMBOLS  PLOTTED  SEE  COMMENTS  AT  START  OF  MAIN  PROGRAM 

C  OR  TABLE  F-8  IN  PLOTIO  IGL  USER  MANUAL  4010C01  (PART  NO.  070-2685-02) 

C  TOR  THE  DEFAULT  ADE  CHARACTERS 

C 

1 50  CONTINUE 

MARK  =  ITYPE-1 
DO  200  IP  =  1 , NPT 

IF(LB(IP) )  CALL  MARKER (X(IP),Y(IP), MARK ) 

200  CONTINUE 

GO  TO  550 
C 

C  PLOTS  SEMI-HISTOGRAM  ( ICHAR=1 7 )  OR  FULL  HISTOGRAM  (ICHAR=18) 

250  CONTINUE 

DO  300  IP  =  1 , NPT 

IF( .NOT. LB ( IP ) )  GO  TO  300 

II =MAX( IP-1 ,1 ) 

I2=MAX( IP/ 2 ) 

I3=MIN( IP, NPT- 1 ) 

I4=MIN ( IP+1 , NPT) 

XDEL1  =  ( X( 12 ) -X( II ) )/2.0 
XDEL2=(X(I4)-X(I3) )/2.0 
X1=MAX(X( IP ) -XDEL1 ,XMIN) 

X  2=MIN ( X ( IP ) +XDEL2 , XMAX ) 

IF( IP.LT.NPT) Y2=Y( IP+1 ) 

IF ( IP. EQ. NPT) Y2=Y ( IP ) 

IF( ICHAR.EQ. 1 7 )  THEN 
CALL  MOVE ( XI , Y( IP) ) 

CALL  DRAW ( X2, Y ( IP ) ) 

CALL  DRAW ( X2 , Y2 ) 

ELSE  IF( ICHAR.EQ. 18)  THEN 
CALL  MOVE( XI , YMIN ) 

CALL  DRAW (XI ,Y(IP) ) 

CALL  DRAW ( X2 , Y ( IP ) ) 

CALL  DRAW (X2, YMIN) 

CALL  MOVE( X2 , Y( IP) ) 

CALL  DRAW(X2,Y2) 

END  IF 

300  CONTINUE 

GO  TO  650 
C 

C  DRAW  LINES 

C  ICHAR  =  20  -  CONTINUOUS  LINE 

C  ICHAR  >  20  -  DASHED  LINE 

C 

350  CONTINUE 

LINE  =  ITYPE-20 
CALL  DASHPT( LINE) 

FIND  FIRST  POINT  WITHIN  FRAME 
DO  400  IP  =  1 , NPT 
IF( LB( IP ) )  THEN 

CALL  MOVE(X( IP) , Y( IP) ) 

I FIRST  =  IP 
GO  TO  425 
END  IF 
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CONTINUE 


400 

C 

C  NOW  PLOT  THE  LINE 

425  CONTINUE 

IF( IFIRST.GE.NPT)  GO  TO  650 

DO  450  IP  =  IFIRST+1 , NPT 

IF(LB( IP ) )  CALL  DRAW ( X  ( IP ) ,  Y  ( IP  ) ) 

450  CONTINUE 

C 

C  RESET  DEFAULT  LINE  TYPE  TO  CONTINUOUS 

CALL  DASHPT(O) 

GO  TO  650 
C 
C 

C  ERROR  BARS 

C 

550  CONTINUE 

DO  600  IP  =  1 , NPT 

IF ( . NOT. LB ( IP ) )  GO  TO  600 

IF(XE(IP) .GT.0.0)  THEN 

XI =MAX ( X ( IP) — XE( IP )  , XMIN) 

X2=MIN ( X ( IP ) +XE ( IP ) , XMAX ) 

CALL  TICK (XI  , Y( IP) ,X2, Y( IP) ) 

END  IF 

IF(YE( IP) .GT.0.0)  THEN 

IF(LOGPLOT. EQ.. FALSE. )  THEN 
Y1 =MAX( Y( IP )-YE( IP ) , YMIN) 

Y2=MIN(Y(IP)+YE(IP) , YMAX) 

CALL  TICK ( X ( IP ) , Y1  ,  X  ( IP ) ,  Y2 ) 

ELSE  IF(LOGPLOT.EQ. .TRUE. )  THEN 

Y1  =  MAX(Y(IP)-LOG10(YE(IP) ),YMIN) 

Y2  =  MIN( Y( IP )+LOG1 0( YE( IP) ) , YMAX) 

CALL  TICK ( X ( IP ) , Y1  ,  X  ( IP ) ,  Y2  ) 

END  IF 
END  IF 

600  CONTINUE 

C 

650  CONTINUE 

RETURN 
END 

C* ********************************************************************* 

C  SUBROUTINE  LOGTRAN 

C********************************************************************** 

SUBROUTINE  LOGTRAN ( X,Y, XE, YE, LY, LYE, NPT) 

C 

C  CONVERTS  DATA  INTO  LOG) 0  REPRESENTATION 

C 

C  INPUT  :  Y, YE  VALUES 

C  OUTPUT  :  LY  -  THE  LOGIO  VALUES  OF  THE  Y  VALUES 

C  :  LYE  -  RETURNED  AS  THE  YE( )  VALUES 

C 

REAL  X( * ) , Y( * ) , XE( * ) , YE( * ) ,  LY ( * ) , LYE( * ) 

INTEGER  NPT 
C 

DO  1000  IPT  *  1 , NPT 
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IF(Y(IPT)  .GT.  0.0)  THEN 
LY(IPT)  =  LOG10(Y(IPT) ) 

ELSE 

LY(IPT)  =  -10.0 
END  IF 
C 

LYE(IPT)  =  YE(IPT) 

C 

1 000  CONTINUE 
RETURN 
END 

C****** ********************************************** ****************** 

C  SUBROUTINE  LINSCALE 

C* ********************************************************************* 

SUBROUTINE  LINSCALE( XMIN , XMAX , XMV , XD IV , IPX , NTX ) 

C 

C  INPUT  :  XMIN,  XMAX  -  THE  MIN,  MAX  VALUES  TO  SCALE 

C  OUTPUT  :  XMV  -  THE  SCALED  MINIMUM  VALUE  (NO.  FROM  0  TO  10) 

C  XDIV  -  THE  SCALED  INCREMENT  BETWEEN  LARGE  TICK  MARKS 

C  IPX  -  THE  POWER  OF  10  MULTIPYING  THE  VALUES 

C  NTX  -  THE  NUMBER  OF  INCREMENTS  ON  THE  AXIS 

C 

C  FIRST  RESET  XMIN,  XMAX  IF  XMIN  GE  XMAX 

IF  (XMIN  .GE.  XMAX)  THEN 

IF  (XMIN  .EQ.  0.0)  XMAX  =  1.0 

IF  (XMIN  .NE.  0.0)  XMAX  =  2.0  *  ABS(XMIN) 

END  IF 
C 

DIFF  =  (XMAX-XMIN) /FLOAT (NTX) 

CALL  XPOW( XMIN, B1 , IP1 ) 

CALL  XPOW ( XMAX, B2 , IP2 ) 

IPX=MAX ( IP1 , IP2 ) 

IF(IP1  .EQ.  0)IPX  =  IP2 
XMV  =  XMIN/1 0. **IPX 
XDIV  =  DIFF/1 0. **IPX 
RETURN 
END 

C* ********************************************************************* 


C  SUBROUTINE  LOGSCALE 

C* ********************************************************************* 

SUBROUTINE  LOGSCALE( YMIN , YMAX , YMV, NTY , IP Y , INC , ISPEC ) 


C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 


SETS  VERTICAL  SCALE  FOR  LOG  PLOTTING 

INPUT  :  YMIN, YMAX  -  THE  MIN,  MAX  VALUES  TO  SCALE 

NOTE  THESE  VALUES  MAY  BE  RESET  BY  LOGSCALE 
OUTPUT  :  YMV  -  THE  SCALED  MIN.  VALUE  (NO.  FROM  1  TO  10) 

NTY  -  THE  NUMBER  OF  INCREMENTS  ON  THE  AXIS 
IPY  -  THE  MIN.  POWER  OF  10 

INC  -  THE  EXPONENT  INCREMENT  FOR  THE  LARGE  TICK  MARK 
ISPEC  =  0  -  SCALES  IN  POWERS  OF  10  ONLY 
ISPEC  =  1  -  SCALES  1 . ,2. ,5.  X  POWERS  OF  10 


INTEGER  IY(  2 ) 
REAL  Y( 2 ) 


C 
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ISPEOO 
INC=1 
YMV=1  .0 

IF  (YMIN  .LE.  0.0)  YMIN=1 .E-l 0 
IF  (YMAX  .LE.  0.0)  YMAX=1.E10 
IF  (YMAX  .LE.  YMIN)  YMAX=YMIN*1 0.0 
CALL  XP0W( YMIN, Y1 , IY1  ) 

CALL  XPCW ( YMAX, Y2, IY2 ) 

IF(Y2.GT.1 .5)  IY2=IY2+1 
NTY=IY2-IY1 

BRANCH  FOR  DIFFERENT  VALUES  OF  NTY 
IF(NTY  .LE.  5)  GO  TO  100 
IF( NTY  .GT.  10)  THEN 
INC-(NTY+9)/1 0 
NTY=( NTY+INC-1 )/INC 
END  IF 
IPY  *  IY1 
YMIN=10.**IPY 
YMAX=10.**(NTY*INC+IPY) 

RETURN 

ISPEC=1  CASE  -  SCALES  1.,2.,5.  X  POWERS  OF  10 
1 00  ISPEC=1 

CALL  XPOW ( YMIN, Y  ( 1  )  , IY ( 1 ) ) 

CALL  XPOW ( YMAX ,Y(2),IY(2) ) 

DO  200  1=1,2 

IF(Y(I) .GT.7 . ) IY( 1) =IY( I)+1 
IF( Y( I) .LE. 1 .5.0R.Y(I). GT . 7 . )  Y(I)=1 .0 

IF( Y ( I ) .LE. 7 . .AND • Y( I ) .GE. 3.5)  Y(I)=5.0 
IF(Y( I) .LT.3.5.AND.Y(I).GT.1 ,5)Y(I)=2.0 
200  CONTINUE 

C 

YMV=Y(1 ) 

IPY=IY ( 1 ) 

YMIN=Y( 1 )*1 0. **IY(  1  ) 

YMAX=Y(2)*10.**IY(2) 

IT1  =0 
IT2=0 

IF( Y( 1 ) .NE. 1.0)  IT1 =( +Y{ 1)+1.1)/3.0 
IF( Y( 2) .NE.1 .0)  IT2=(+Y(2)+1 .1 )/3.0 
NTY=3*(IY(2)-IY(1  ) )  +  (IT2-IT1) 

IF( Y( 1 ) . EQ . 2 • • AND  «Y(2) . EQ . 5 . . AND • I Y ( 1 ).EQ.IY(2) )  NTY=1 
IF( Y( 1 ) . EQ • 5 . • AND . Y ( 2 ) . EQ • 1 ,.AND.IY(1 ).EQ.IY(2)-1 )NTY=1 
RETURN 
END 

Q* ******* ************************************************************** 

C  SUBROUTINE  XPOW 

C* ******* A**************** ************************************ ********* 

SUBROUTINE  XPOW(XVAL,XNORM, I EXP) 

C 

C  INPUT  :  XVAL  -  ANY  REAL  NO. 

C  OUTPUT  :  XNORM  -  10  RAISED  TO  THE  MANTISSA  OF  LOGIO(XVAL) 

C  I EXP  -  THE  CHARACTERISTIC  OF  LOGIO(XVAL) 

C  NOTE  DOUBLE  PRECISION  USED  TO  AVOID  ROUNDING  ERRORS 
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C 

DOUBLE  PRECISION  DA,  DB,  DG 
C 

DA  =  DABS (DBLE(XVAL) ) 

IF(XVAL  .EQ.  0.0)  THEN 
XNORM  =  0. 

I EXP  =  0 
RETURN 
END  IF 
C 

DG=DLOG1 0(DA)  +  1.D-6 
I EXP  *  IDINT(DG) 

IF( DG  .LT.  0.D0)  IEXP  =  IEXP-1 
DB  =  DA  /  DBLE( 1 0.0** IEXP) 

IF( DB  .LT.  1 .DO)  DB  =  1.D0 
IFCXVAL  .GE.  0.0)  XNORM  =  SNGL(DB) 

IF(XVAL  .LT.  0.0)  XNORM  =  -1.0  *  SNGL(DB) 

RETURN 

END 

q* ********************************************************************* 

SUBROUTINE  TICK 

DRAWS  TICK  MARKS  OR  LINE  SEGMENTS  BETWEEN  TWO  SPECIFIED  POINTS 
*************************************************************  ********* 
SUBROUTINE  TICK  ( XI , Y1 ,X2, Y2 ) 

REAL  XI ,Y1 , X2, Y2 

XI, Y1  ARE  COORDS.  OF  START  PT,  X2,Y2  COORDS  OF  STOP  PT. 

CALL  MOVE (XI ,Y1 ) 

CALL  DRAW(X2,Y2) 

RETURN 
END 
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